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Abstract 

We utilize the DVCS asymmetry measurements of the Hermes collaboration for access to Comp- 
ton form factors in the deeply virtual regime and to generalized parton distributions. In par- 
ticular, the (almost) complete measurement of DVCS observables allows us to map various 
asymmetries into the space of Compton form factors, where we still rely in this analysis on dom- 
inance of twist-two associated Compton form factors. We compare this one-to-one map with 
local Compton form factor fits and a model dependent global fit. 
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1 Introduction 



As is well known, at the beginning of the last century Compton scattering played a crucial role 
in the debate on the dual nature of light. Namely, if light scatters off an electron, energy is 
transferred to the electron and so the wavelength of the scattered light increases. In addition 
to this Compton effect revealing the particle aspect of light [T], the Compton scattering process 
has now many applications in material science, medicine, astro- and particle physics, where in 
particular, it is utilized to probe the electron wave function in matter. In hadronic physics 
the Compton scattering process is used to reveal static properties of the nucleon in terms of 
electric and magnetic polarizabilities. Its generalization into virtual Compton scattering [2] offers 
a supplementing window to electromagnetic form factor studies, which have been intensively 
performed in the last six decades. At larger photon virtualities > IGeV^, the Compton 
scattering process probes the partonic content of the nucleon and gives access to the so-called 
generalized parton distributions (GPDs) [31 HI [S]. In particular, it was suggested to measure 
by means of deeply virtual Compton scattering (DVCS) the quark orbital angular momentum 
[6] and to access the transverse distributions of partons [71 [8]. With these goals in mind, many 
experimental and theoretical activities were started in the field of hard exclusive processes, 
where we consider the measurements of (DV)CS observables as important as the measurements 
of electromagnetic form factors. 

Measurements of DVCS have taken place at Hera and Jefferson Lab since the turn of the 
millennium, with the aim of understanding the decomposition of the nucleon spin. The first 
observed DVCS candidate events at high energy were announced by the Zeus collaboration 
using a collider experiment at Hera in 2000 [9] and few months later for fixed target kinematics 
[To] . The first DVCS measurements in fixed target experiments were published simultaneously 
by the Hermes and Clas collaborations a year later [Til [12] • Both experiments measured the 
beam spin asymmetry arising from the interference term between DVCS and elastic scattering 
with bremsstrahlung in the scattering amplitude. The latter process, called the Bethe-Heitler 
(BH) process, has the same initial and final states as DVCS {ep — )■ ep^) and at fixed target 
kinematics typically the interference term in the scattering amplitude is more experimentally 
accessible than the pure DVCS term [THj. Both Hermes and Clas utilized the fact that the 
large BH contribution in the interference term at ~ 2 GeV^ and —t ~ 0.1 GeV^ amplifies the 
contribution from the more interesting DVCS process. 

Since then, large experimental effort has been expended to measure various observables 
in the electroproduction of photons and mesons. The phenomenological challenge is now to 
describe these data in terms of GPDs, which requires a flexible model. Thereby, we meet a more 
elementary problem, in particular for DVCS: the number of Compton form factors (or helicity 
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amplitudes) is usually larger than the number of observables at a given kinematical point. One 
must therefore rely on model assumptions or hypotheses which means that, independently of 
the applied method or framework, a theoretical bias cannot be avoided in analyzing the present 
available world data set. This may even lead to a qualitative misinterpretation of the data, 
where the apparent influence of any particular GPD model set may be determined more by the 
initial theoretical approach to the problem than by the observed data. 

In this paper we will concentrate on the results from the Hermes collaboration, which had 
both electron and positron beams available and is currently the experiment that has delivered 
the most complete set of DVCS observables. In Sect. [2] we introduce the DVCS observables 
and we give details on the Hermes experiment. In Sect. Owe consider the extraction of CFFs 
at given kinematical points from the Hermes measurements as a map of random variables. 
Additionally, we utilise the regression approach and use the Hermes data to access CFFs by 
least squares fitting. We also present a global GPD model fit that additionally includes Hera 
collider and Jefferson Lab measurements and conclude on what we have learned. Finally, we 
summarize and provide an outlook. 

2 HERMES measurements of DVCS observables 

2.1 Definitions of photon leptoproduction observables 

The DVCS process enters as a subprocess in deeply virtual photon leptoproduction and its 
amplitude can be decomposed in terms of twelve independent helicity amplitudes given in some 
reference frame. Alternatively, we might parameterize the DVCS amplitude in Lorentz-invariant 
Compton form factors (CFFs), which are defined in analogy to GPD definitions and are called 

I'H, "H, , Jeff G I'Heff, i^eff, "Heff, i^effj , J^T G jT^T, ^^T, "Ht, , (1) 

where we separated them into twist-two related CFFs J-", twist-three related CFFs J-'cs, and 
transversity CFFs J-t. Unfortunately, different conventions, which differ by power suppressed 
contributions, are used in the literature. As noted previously, the DVCS subprocess is accom- 
panied by the Bethe-Heitler bremsstrahlung process, the amplitude of which is, at leading order 
in the QED fine structure constant, real- valued. Moreover, it is entirely parameterized in terms 
of the electromagnetic form factors, which for the nucleon are well-known for the kinematics 
of interest. The interference term of both subprocesses is charge-odd and depends linearly on 
the CFFs, while, the BH-squared, depending only on electromagnetic form factors, and DVCS- 
squared, given as bilinear form of CFFs, amplitudes are charge-even. The differential cross 
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Figure 1: Reference frame adopted by the Hermes collaboration [T^. 



section, for the most general setup, is five-fold and we write it in terms of the Bjorken scaling 
variable xb, the negative virtuality squared of the intermediate photon = —q"^, the mo- 
mentum transfer square in the t-channel, and two relative azimuthal angles 0, between lepton 
scattering plane and reaction plane, and (p, between lepton scattering plane and transverse spin 
polarization vector, as in [Ti] : 

(2) 

where a = e^/47r is the electromagnetic fine structure constant and Mp is the proton mass. Here, 
and in the following, we neglect the electron mass and define the azimuthal angle in an reference 
frame where the photon travels along the 2-axis and the x-component of the incoming electron 
is positive, see Fig. [T] This frame has been adopted by the Hermes collaboration. Compared 
to the BMK conventions in [T3], in which the photon travels in opposite direction of the 2;-axis, 
for the frame used by Hermes we have the relations 

0BMK ^ ^BMK _ ^BMK ^ 0, - 0e (3) 

^BMK ^ ^BMK _ ^B..K ^ ^ _ ^ ^ ^ = 0,-08, (4) 

where (pN, (p-y = (pN + tt, 0e = 0, and $ = 0s are the azimuthal angles of the nucleon, photon, 
electron, and the transverse polarization vector, respectively. 



BHI 



(xb, Q'^,t,(j),(p) 
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The Fourier expansion of the various amphtude squares may be written in analogy to the 
notation of [13] as 
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'^bhI'= ,2 r,A cr + J2^n''cos{n(f)) + rs^''sm<f)\ , (5) 
—tx^w[(p) I ^-^ I 



±e 



n=l 
3 



" _^^2 ,„u^ i ^0 + K cos(n0) + si sin(n0)] ^ , (6) 




Here, the + (— ) overall sign of the interference term ([6]) refers to an electron (positron) beam, 
l/w{(f)) oc 1/7^1(0)^2(0) is the product of scaled BH propagators given in [2], which is expanded 
in the even harmonics of up the second order, 

w{(f)) = 1 + Wi cos(0) + W2 cos(20) . (8) 

Note that we included a generic kinematical overall factor \/ —tjiy^Q?^ from the interference 
term ([6]) in the definition of the Fourier coefficients. Moreover, if the differential cross section 
([2]) is weighted with w(0), its Fourier expansion w.r.t. the azimuthal angles becomes finite. 

Having both electrons and positrons with which to make measurements, one can access the 
interference and DVCS-squared terms by means of cross section differences and sums 

LOqed 1 1 

Measuring such charge-odd cross sections with different lepton and nucleon polarizations, one 
can, in principle, reach an almost complete decomposition in terms of the twelve CFFs ([1]). 
Namely, in the deeply virtual regime, the first and second order harmonics of the weighted 
differential cross section differences are dominated by linear combinations of twist-two and twist- 
three related CFFs J-" and TcS respectively, while the third harmonic is related to the transversity 
CFFs J-T. In the first two cases the number of possible observables allows access to both the 
imaginary and real parts of the corresponding form factors, while in the latter case only one 
combination of iHe J-r and three combinations of JmJ-r can be accessed. Moreover, in this 
approximation the constant terms are governed only by twist-two associated CFFs J-" and, hence, 
such measurements allow for an experimental consistency check of the underlying formalism, 
given in some approximation. 

The charge-even combination may serve also for an experimental consistency check or may 
be used to access the desired CFFs in a different manner. Most importantly, in single spin-flip 
experiments, the BH term drops out in the considered order of QED and so such measurements 
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give direct access to bilinear CFF combinations. In double spin-flip experiments one needs to 
subtract the BH-squared term, which may be feasible at small xb where an effective 'pomeron' 
behavior may overwhelm the contributions of the BH subprocess. 

Thus far, only Hera experiments at Desy had both charges of lepton beams at hand. The 
Hermes experiment achieved measurements with polarized electron and positron beams for 
longitudinal and transversal target spin polarizations and thus the following asymmetries where 
extracted. 



Single beam-spin asymmetries in the charge-odd and charge-even sectors: 

[da^ - daX] - [da:! - da!:] 



^LU,l(0) 
^LU,DVCs(0) 



da^ + daX + da^ + dat 
[da']^ - daX] + [da:! - da^] 



da':! + da\ 



da_ 



dal 



(10) 
(11) 



where — )■ (4—) denotes electron (— ) or positron (+) polarization along (opposite to) the 
beam direction. 



Beam-charge asymmetry: 



da^ — da^ 
da^ + da. 



(12) 



Single and double longitudinal target-spin asymmetries with a positron beam: 

. X daT — da^ , , 

AuL,+(0) = , , ; , (13) 

dax + da^ 

[daX^ + da^^\ - [da^^ + daX^\ 



da\ 



da. 



da. 



da*. 



Note that, contrary to Hermes notation, <^= (^) denotes proton polarization opposite to 
(along) the positron beam direction. 



Single transverse target-spin asymmetries in the charge-odd and charge-even sectors: 



^UT,l(0, V^) 
^UT,DVCS(0, V^) 



dat. — da. 



dal - da 



dal + dal + rfa'^ + da 



da^ — da 



da^ — da 



dal + dal + dal + da 



(15) 
(16) 
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Double longitudinal beam and transverse target-spin asymmetries in the charge-odd and 
charge-even sectors: 



^LT,i(0,V5) = — T-^r7 

n.cT . + 



- da^^ + daX^ 



(17) 



■LT.evenl 



rf^xr^ + rfcxr^ + da^^ + da^"^ - da^^ + da^^ 



(1^ 



All these observables can be theoretically expressed in terms of weighted cross sections. However, 
since the denominator given by the unpolarized cross section depends on the azimuthal angle 
0, they are given by a series of harmonics rather than finite sums. Nevertheless, in the 
kinematics where the BH process dominates, the zeroth and first harmonics of these asymmetries 
are roughly determined by the linear combinations of twist-two associated CFFs J-" that enter 
the interference term. Note also that, except for the longitudinal target spin asymmetries f|T3l) 
and f|T^ . the denominators are expressed by the charge-even cross section. Utilizing the charge 
asymmetry f[T^ . we find in general that an asymmetry measured with a positron (or electron) 
beam can be expressed by charge-odd (interference term) and charge-even (a possible squared 
BH term plus a DVCS-squared term) expression: 

A..,i(0) , A..,bh(0) + A.. ,DVCS (0) 



A..,+ (0) 



(19) 



1 + Ac(0) 1 + Ac(0) 

Assuming that the BH amplitude overwhelms the DVCS one, we may drop the DVCS induced 
part and we approximately have 



A..,i(0) ^ A..,+ (0) [1 + Ac(0)] - A..,bh(0) 



(20) 



where the BH-associated asymmetry A... bh(0) drops out in single-spin asymmetries. We add 
that the double-spin asymmetry in the charge-odd sector is defined in terms of cross sections as 

[da^^ + da^^] - [da^^ + da^^] - [da^^ + da^^] + [da^^ + da^^] 



da^^ + da A 



da' 



+ da^^ + da_ 



da' 



da^^ + da: 



(21) 



Let us finally remind the reader here that asymmetries are expanded in an infinite Fourier 
series. Hence, the method for the extraction of CFFs, outlined above for cross section differences, 
has to be modified. For illustration let us consider here the lepton beam-spin (fTO|) and beam- 
charge (fT2|) asymmetries that offer access to the imaginary and real parts of CFFs respectively. 
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Substituting the cross section ([2]) and the harmonic expansions (15|6|7|) into ( ITOj) and ( !T2|) . we find 
the azimuthal angle dependencies of the charge-odd electron beam spin and charge asymmetry 
respectively: 

. (,^ sin(0) + gf sin(20) 

^LU,il0j - ^ 7^ — , IT^TT?^^ l^^^ 



E:=o cJ?^ cos(n0) + -^^(0) E:=o c^^^' cos 



^ . ^ eg + cos((/>) + cf cos(20) + cf cos(30) 



E:=o cJ?^ cos(n0) + ^^(0) E:=o 2^^^' cos( 



where the odd and even Fourier coefficients and of the interference term are linear functions 
of the imaginary and real parts of the CFFs and the even DVCS Fourier coefficients cj^^*"^ 
are bilinear in all CFFs. In the 1/Q expansion, the first harmonics of the interference term 
dominate and are governed by the CFFs J-", the second harmonics are relatively suppressed 



by (tmin — t)/Q'^ and determined by the CFF J-'^g while the third harmonic of the charge 
asymmetry fl25]l is given by the CFFs J-t and, to leading-twist accuracy, is predicted by the 
gluon transversity GPD. In reality, however, all harmonics of the interference term are functions 
of all twelve CFFs [T6| [T7] . Note that the third odd harmonic is absent and hence the imaginary 
part of the transversity CFF combination cannot be accessed from the interference term. In 
principle, this missing information is contained in the DVCS-squared term. Besides the mixing 
of the various CFFs for a given harmonic, the projection 

^sinM) ^ 1 sin(n0)ALu,i(0) , (24) 

^cosin^p) ^ 1 cos{n^)Ac{4>) for n > 0, A'°'^''^^ = — /"rf</.ylc(0) , (25) 

yields an additional contamination due to the dependence in the denominator. 



2.2 Experimental details 

Exclusive photon events at Hermes were selected if having exactly one photon and one lepton 
track detected within the acceptance of the spectrometer. The event selection was subject to 
the kinematic constraints 

1 GeV^ < < 10 GeV=^ , 0.03 < xb < 0.35 , -t < 0.7 GeV^ , W"^ > 9 GeV^ , and z/ < 22 GeV , 

where W is the invariant mass of the 7*p system and u is the energy of the virtual photon in 
the target rest frame. The polar angle between the directions of the virtual and real photons 
was required to be within the limits 5 mrad < 6^*^ < 45 mrad. 

An event sample was selected requiring that the squared missing- mass = {q + Mp — g')^ 
of the — )■ e'7X measurement corresponded to the square of the proton mass, Mp, within the 
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1 • 

bin no. 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


-t [GeV'] 


0.03 


0.1 


0.2 


0.42 


0.1 


0.1 


0.13 


0.2 


0.08 


0.1 


0.13 


0.19 




0.08 


0.1 


0.11 


0.12 


0.05 


0.08 


0.12 


0.2 


0.06 


0.08 


0.11 


0.17 


[GeV^] 


1.9 


2.5 


2.9 


3.5 


1.5 


2.2 


3.1 


5.0 


1.2 


1.9 


2.8 


4.9 



Table 1: Kinematical values of three times four Hermes bins from ref. [13], labeled as #1, • • • , i^l2. 

limits of the energy resolution of the Hermes spectrometer (mainly the calorimeter). Recall 
that q is the four-momentum of the virtual photon, p is the initial four-momentum of the target 
proton and q' is the four-momentum of the produced photon. The "exclusive region" was defined 
as -(1.5 GeV)^ < M| < (1.7 GeV)^. This exclusive region was shifted by up to 0.17GeV^ for 
certain subsets of the data in order to reflect observed differences in the distributions of the 
electron and positron data samples. A systematic uncertainty contribution is assigned for this 
effect. This event sample selection technique was used as it allows for the most complete set of 
DVCS observable measurements to be considered; although results using the measurement of the 
recoiling proton have been released by Hermes [18] , this technique is only currently available for 
beam-spin DVCS measurements and was never utilised with a polarised target. Other systematic 
uncertainty contributions arise from potential misalignment of the spectrometer, acceptance and 
smearing effects, and the inclusion in the data set from misidentified semi-inclusive deep inelastic 
scattering events. The latter is mitigated somewhat using a correction procedure. The sizes of 
the systematic uncertainties are estimated using Monte Carlo techniques. 

The data sample in the exclusive region contains events not only involving the production of 
real photons in which the proton remains intact, but also events involving the excitation of the 
target proton to a A"^ resonant state ("associated production"). This is a consequence of using 
a missing-mass event sample selection technique as noted above; the calorimeter resolution for 
measuring the produced photon does not allow separation of the resonant events from the rest of 
the data sample. No systematic uncertainty is assigned for the contribution from these events; 
they are treated as part of the signal. A Monte Carlo calculation based on the parameterisation 
from ref. [IH] is used to estimate the fractional contribution to the event sample from resonant 
production in each kinematic bin; the uncertainty on this estimate cannot be adequately quan- 
tified because no sufficiently precise measurements have been made in the Hermes kinematic 
region. The results of the estimate, called the associated fractions and labelled "Assoc. frac- 
tion", are shown in refs. [131 [151 [201 [211 [22]- The method used to perform this estimation is 
described in detail in ref. |13] . 

The Hermes measurements were presented in terms of four or six bins times three one 
dimensional projection in —t, xb, and (see Tab.[T|). The projections in the kinematic variables 
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DVCS asymmetries 


"^'^^^^ 






Coniptonform factors 


• data ^ CFFs 






DVCS asymmetries 







HERMES 


— / = 0.12 GeV~ 




[BMIO/BMKOI convenlions) 


xs = 0.097 

— / — U. 1 Z Kjc V 




HERMES 


A'a = 0.097 
—1 = 0.12 GeV~ 






Q- = 2.51 GeV^ 






= 2.51 GeV^ 






= 2.51 GeV^ 






o CFFs => data 












o CFFs ^ data 


sin (ip) 


i-o-i 




Im 'H 












sin U] cos M 






Im & 






sin (^tp'j cos (O (hj 
'^UT,DVCS 














1 -o 




-COS {tp) sin (<p) 






Im £ 






cos (ip) COS (O S) 
'^LT,BH+DVCS 






1— ( 
1- 










) 


-COS M 




• 

O 


■Re 




H 


COS (00) 


m 

G 
















cos (0 ip) 






1- -o- - 














COS (v) COS W 






2 






cos [tp) COS (0 1^) 






K - e - 




1 • 








-OH 


sin (v) sin (0) 






-xKe S 

2 






sin (^) cos (00) 









© - H 




l_l 1 1 1 1 1 1 1 1 1 1 










' ' Ij I I I I I I I I I I I I I I I I I I I I J ' ' 

-0.3 -0.2-0.1 0.0 0.1 0.2 0.3 -20 -10 10 20 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 

asymmetry value Compton form factor value asymmetry value 



Figure 2: Hermes measurements (solid circles) of interference dominated asymmetries in the twist- two 
sector (left) and BH/DVCS-squared dominated as well as charge odd twist-tree asymmetries (right) at 
the overall mean values XB = 0.097, t = -0.119 GeV^ and = 2.51 GeV^ In the middle panel the 
resulting CFFs from a linearized map of the interference dominated asymmetries are displayed. The 
empty circles, shown in the left and right panels, arise from the inverse linearized map of extracted 
CFFs back to asymmetries. 

are each correlated; the and x-q projections are very highly correlated as a consequence of 
the experimental design. For our analyses we employ the four bin data, where the four-binned 
version of the six-binned data for the beam spin and charge asymmetries from the Hermes 
measurements, published in ref. [22], were provided to us by the Hermes collaboration. We 
neglect the small differences in the kinematical means of different asymmetry measurements, 
which are very much within the experimental uncertainty and take the values that are listed in 
Tab. [H labeled as bin 7^1, ■ ■ ■ , #12. Alternatively, we take the kinematical values that are given 
in the publications, assuming that the asymmetry values do not vary within the small changes 
in the mean kinematic values. The selected set of fourteen observables out of the thirty-four 
asymmetry measurements, which we will use, are shown in Fig. |2] as solid circles at the overall 
mean values = 0.097, t = -0.119 GeV^ and = 2.51 GeV^ 

3 Analysis of HERMES data 

In Sect. 13.11 we first illustrate which information one may extract from the Hermes measure- 
ments, cf. Fig. 121 if one assumes that effective twist-three and transverse photon helicity flip 
effects are absent. Employing the BMK formulae set [H], we will also show that data are not in 
contradiction with this hypothesis. The outline of the remainder is as following: in Sect. 13.21 we 
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give a review of the methods used to analyze DVCS data and introduce the (non-)hnear mapping 
method. In Sect. 13.3.11 this method is then apphed to the Hermes data and in Sect. 13.3.21 we 
analyze the same data with local CFF fits. Finally, in Sect. 13.41 we present a global CFF fit 
with a simple GPD model used previously in global fits to DVCS data off unpolarized protons. 
Finally, we discuss the implications of our analyses for GPD model building. 



3.1 Asymmetries in terms of CFFs 

Relying on the dominance of the BH-squared term over the DVCS-squared term and approximat- 
ing the BH-squared term by its dominant constant harmonic, the following linear combinations 
of CFFs enter in the various asymmetries (here we set the minimal value of —t to zero): 

• The first harmonics of single beam spin ( fTOj) and beam charge (fT2|) asymmetries. 



are approximately given by the imaginary and real part of the linear combination 

CUi^) ^ F,n - -^F,S + Y(i^i + F2)n . (28) 

The first harmonics of single longitudinal target spin asymmetry f|T3l) and double longitu- 
dinal target spin asymmetry ([1] 



^ul!+^ oc -^^y-^xa;B3mKp(^) + -.-]+DVCS2-term, (29) 

^LL^*^ oc - y^^^ ^EL X Kp(^) + ■■■]+ BH^-term + DVCS^-term , 

(30) 

are approximately governed by the imaginary and real part of the linear combination 

C[p(^) ^ F,1l - i^F, + ^F,) ^(Fi + F,)n . (31) 

The first even harmonic of the single transverse-target spin asymmetry f|T5|) and first odd 
harmonic of the double transverse-target spin asymmetry (ITTl) in the charge-odd sector. 



^si„(.)cos(i,) ^ -y^J-^ X XB 3m [4p„(-F) + ■■■], (32) 



A^^-^^ oc ^ X [4p„(-F) + ■■■], (33) 
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are approximately dominated by the imaginary and real part of the hnear combination 



4p_(^) 



-t 



4M2 L 



2F2H - 2F,S + xb(Fi + F2)y^ 



(34) 



The first odd harmonic of the single transverse-target spin asymmetry f|T5|) and first even 
harmonic of the double transverse-target spin asymmetry f|T7j) in the charge-odd sector, 



^cos((/3) sin(l(/)) 



1/ 



M2 

p 



2-y \y^Q^ 



X 



XBam[4p+(J^) + ---] 



(35) 



.cos((p)cos(l<}!)) 
/Ilt T OC -1- 



M2 

p 



2 -2?/ + 2/2 W ^232 



X 



XB [4p+(^) + ---] 



(36) 



are approximately dominated by the imaginary and real part of the linear combination 

-t 



4M2 

p 



2F2H 



+ ■ ■ 



(37) 



As one realizes from the relations among observables and CFFs listed above, there are six 
linear combinations of CFFs: three for the imaginary parts fl26ll32ll35p and three for the real 
parts fl27ll33ll36p . Unfortunately, the longitudinally polarized target double spin flip experiment 
has been performed only with positron beam and, hence, the asymmetries (l29ll and (l30l) are 
contaminated by twist-three contributions from the DVCS-squared term. Moreover, the latter 
term also depends on the first harmonic of the BH-squared term. Since the DVCS-squared 
contributions are expected to be relatively small and the first BH harmonic is suppressed for 
Hermes kinematics, the single spin and double spin flip asymmetries may be used to access 
the imaginary and real parts of the CFF combination (13T|) . To get rid of these contaminations, 
we would need longitudinal single and double spin flip measurements in the charge odd sector 
which were unfortunately not performecj^. 

Further twist-two dominated observables are the single transverse proton spin asymmetry 
and the double longitudinal-transverse spin asymmetry in the charge even sector: 

• The zeroth harmonic of the single transverse target spin asymmetry, which is governed by 
the bilinear combination 

1 - y i-tf/' 



A 



sin((p) cos(0</<) 
UT,DVCS 



OC + 



X x| 3m 



(38) 



4(2-2/) Mp2/2g2 

The zeroth harmonic of the double longitudinal-beam and transverse target spin asymme- 
try, which is dominated by the bilinear combination 

2/(1-2/) (-t)3/2 



,cos(i^)cos(0()!>) 
^LT,BH+DVCS ^ "T 



4(2 -22/ + 2/2) 2/' Q 



2 X x| 9ie 



Xb 



(39) 



-'^For the history of the Hermes experiment see http:/ /www-hermes.desy.de/hedt/seminar.html. 
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The asymmetry in ( 138|) is a small quantity if the CFFs l-i and "H are nearly in phase with £ 
and respectively. Interestingly, in asymmetry f l39p . CFFs with different parity are combined. 
We emphasize that for this longitudinal-transverse double spin asymmetry ( l39l) the BH-squared 
term is, at Hermes, kinematically suppressed: 

^cos(vp)cos(0(^) ^ ^cos(vp) cos(0</)) 

^LT,BH+DVCS ~ ^LT.DVCS • I 

Similarly to the single spin asymmetry ^ut dvcs^'^'' ' ^^^^ double spin asymmetry is therefore 
sensitive to the DVCS-squared term; this is demonstrated by the empty rectangles in the middle 
row of Fig. m that show the pure BH contribution. As mentioned above, the constant terms 
for charge odd asymmetries, which are relatively suppressed by 1/Q w.r.t. the first harmonic, 
may also be dominated by twist-two associated CFFs; however, they may suffer from a larger 
contamination of effective twist-three CFFs, which are neglected here completely. Here we have 
available the real part of two linear combinations and the imaginary part of one: 



(41) 



The zeroth harmonic of the beam charge asymmetry ( fT2l) 



yQ 



2-2y + y' 



The zeroth harmonic of the single and double transverse target flip asymmetries (fTSi) . and 
(HID 



,(42) 
.(43) 



Apart from terms that are kinematically suppressed at small x-q and large y, we realize that these 
asymmetries are already expressed in terms of CFF combinations (!28| [34| [37|) that determine the 
first harmonics. Hence these three asymmetries, measured by HERMES, may roughly expressed 
in terms of the first harmonics: 

^cos(o0) _ ' ^~y jr^<^'t>) C44) 

.sm((p)cos(09i) 2 — y ^ sin(y) cos(10) /.p-x 

^UT,i g-^7f^^UT,i , I45j 

.cos{v5)cos(O0) '\/~^ 2 — y ^cos(y)cos(10) /.^x 

^LT,i g-7f^^UT,i , (46) 

where a fourth relation of the same form exists for the double longitudinal spin asymmetry, see 
(130|1 and (148|1 below. Note that the accuracy of these crude relations can be drastically improved 
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Figure 3: Expectations (squares, sliglitly sliifted to tire left) from tlie approximations (|44H46p and values 
extracted from the one-to-one map (stars, slightly shifted to the right) are compared to the Hermes 
measurements (circles) of ^4^°'^^'^''^) (top), ^"^(^'1') (middle), and j\^(v>) cos(Oct>) ^^Jq^j^^ asymmetries 

in 12 bins by means. 



by taking the higher BH harmonics into consideration. The relations among the charge odd 
asymmetries pi] - [46|) offer an experimental consistency check of the underlying formalism, which 
is shown in Fig. |3l It can be seen that the expectations (squares) of our rough approximations 
|6]) are mostly satisfied, where the largest deviation of ~ 2a appears in the 10**^ bin of 



A 



sm{ip) cos(0(/)) 
UT,I 



. Unfortunately, the double spin asymmetries A- 



cos{ip) cos(0(^) 
LT,I 



and A- 



cos(ip) cos(l(^) 



LT,I 



from limited statistics. 

Another experimental test of the BH dominance is provided by: 

• The zeroth harmonic of the longitudinally double spin flip asymmetry (fT4|) 

jCos(0</i) Acos{0(f>) tcos{0(f>) tcos{0(f>) 

^LL,+ ~ ^LL,BH + ^LL,I + ^LL.DVCS ' 



which is decomposed in its charge odd part 
2-y -t 



,cos{0<f>) 



2-2y + y^ yQ^ 



CUT) + 



suffer 



(47) 
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Figure 4: Hermes measurements of the double spin flip asymmetries (top row) and 

^LTM+DVCS (middle row), as well as of the transverse single proton spin asymmetry Dvcs"'^^ 
(bottom row) in the charge even sector are displayed as circles. The measurements are compared with 
pure BH contributions (empty rectangles, slightly shifted to the left), expectations from the approxi- 
mation (I50p (squares, shifted to the left) and values extracted from the one-to-one map (stars, slightly 
shifted to the right). 



and the charge even part, where the DVCS induced asymmetry part can be safely dropped. 



.cos(0(/.) ^ 
^LL,BH-|-DVCS ~ 



A 



cos(0(/>) 
LL,BH 



(49) 



For Hermes measurements, the relations 



yield the expectation 



^cos(O0) 



cos(0(/)) 
LL,BH 



y 



^C0S(1()!)) _ .cos(10) 



which is visualized in the upper panels of Fig. HI It can be seen that the expectations (squares) 
of our rough approximation f l50|) are satisfied, except in the third and sixth bins. Note that if 
we take these bins literally, we have to conclude that in these two cases the DVCS amplitude 
overwhelms the BH amplitude and so the sign expected from the BH term changes. Perhaps 
more realistically, we may view these two bins simply as a statistical fluctuation on the ~ 2a 
level. 
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3.2 Methods for analyzing DVCS data 

The amplitudes of the harmonics for the various asymmetries, extracted by the Hermes col- 
laboration, can be analyzed in various manners, e.g., one may consider any of the following 
strategies. 

i. CFFs are locally extracted from a complete set of observables by a map. 

ii. CFFs are locally extracted by means of least squares, likelihood, or neural network fits. 

iii. Comparing measurements with model predictions. 

iv. Global model fits based on the least squares method, maximum-likelihood estimation, or 
neural networks. 

The term local refers here to a three dimensional kinematical data point (bin) given by a xb, 
t, and value. The first method has been proposed and discussed in [Tl], however, it has 
not been applied so far. Local extraction of the imaginary and real parts of CFFs (which, for 
brevity, we will call sub-CFFs) by means of least squares fits were applied to Hermes and JLAB 
data [23l|24l[25]. However, although called model- independent, these analysis are biased, too, 
by utilizing the hypothesis that seven sub-CFFs determine all harmonics and satisfy some model 
constraints. Both methods require that the observables have been measured in a common set of 
kinematic bins, where it is assumed that mean values of the kinematical variables are identical 
for all of the measurements. This, however, can be only approximately true. The kinematical 
condition can be relaxed if one employs neural networks. Presently, they have been employed to 
access the imaginary and real part of CFF l-L from Hermes data on unpolarized protons |26j . 

Model predictions for DVCS have existed for over a decade. In the beginning, it was rather 
popular to use Radyuskin's double distribution ansatz (RDDA) for GPDs, given at t = |27] . 
In this method, one uses information on parton distribution functions (PDFs) and form factors 
to build the model, while the GPD arises from a skewing operation that is governed by a profile 
function and the width of this concave function is controlled by a "6" parameter. Such models, 
e.g. given in [28], have been privately distributed as numerical code, referred according to the 
author names of the paper [21] as VGG, and it is often used to compare data with model 
predictions. However, when this is done, authors often neglect to give full details of the code 
version and precise details of the underlying model, e.g. the set of PDFs that are used to constrain 
the GPDs. The same underlying double-distribution ansatz was also used in the BMK model 
[13] , in the numerical code of Freund-McDermott [3U] , and the Goloskokov-KroU model [211 E2] , 
used in an adjusted hand-bag model framework for the description of deeply virtual meson 
production. Certainly there are differences in all these models, however, since these models imply 
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H dominance and H is inherited from a unpolarized PDF parametrization and sometimes from 
the t-dependence of electromagnetic nucleon form factors, it is not very surprising that all these 
models applied to DVCS off unpolarized protons show similar features in a LO analysis. When 
describing to some extent present DVCS measurements with a RDDA based GPD model, one 
may use a small skewness effect, use in a LO analysis NLO PDF parameterizations and neglect 
GPD evolution (important in HI /ZEUS Collider kinematics). Aiming to describe DVCS data, it 
has been also suggested to build "flexible" GPD parameterization by violating their spectral (or 
polynomiality) property [33] (see also [31], where the claim of satisfying polynomiality is simply 
not true). This is an entirely unacceptable recipe (see comments in [35] and [36]) and so the 
phenomenological aspects of such models cannot be discussed. 

Since it was widely realized in the past that most variations of the RDDA usually overestimate 
the size of beam spin asymmetries and do not describe some other observables, it was necessary 
to build flexible GPD models which can be used in numerical fast fitting routines [371 138] , 
where GPDs are modelled in terms of their conformal moments ^H] (related representations 
were suggested in [1011111112] )• Thereby, one tries initially to parameterise the partonic degrees 
of freedom in such a manner that they are accessible from the experimental measurements. We 
will provide more details in Sect. 13.41 

We emphasize that one major impediment to the understanding of DVCS data in terms of one 
specific double-distribution ansatz is the lack of exclusivity in the data. As noted in section 12.21 
there is an estimated average of 12% contamination of the Hermes data by processes involving 
an intermediate A-resonance. While Hermes has published a beam-helicity measurement from 
a very pure sample [TS] that indicates that this contamination only acts as a small dilution 
of the asymmetry magnitude at higher values of —t (e.g. —t > 0.25 GeV^), there is sufficient 
uncertainty that one cannot make any definitive statements on the specific form of a double- 
distribution ansatz. 

In the next section we outline the method of mapping the asymmetries to CFFs. Thereby we 
consider the asymmetries as normally distributed random variables and consider the extraction 
problem of CFFs simply as a map of random variables, where the functional dependence of 
asymmetries on CFFs is provided by the theory. In fact, for normally distributed random 
variables the formalism we employ is very well known from (non-)linear regression and we need 
only to discuss the maps of means and variances rather than probability distributions. 

3.2.1 Maps of normally distributed random variables 

Let us first clarify what we call the twist-two dominance hypothesis and the additional approx- 
imation that arise from a linearization procedure. For instance, the first harmonic of the beam 
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spin asymmetry ( 122|) in the charge odd sector reads 
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sin(l(^) 
LU,I 
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(i</) sin(</)) 



(51) 



s\ sin(0) + X S2 sin(20) 



E:=o cos(n0) + p^^(0)cDVCS + X ^^(0) ELi cos(n<^ 



'He 



=Wt=0 



where x ■ ■ ■ indicate the expressions that are induced by twist-three and transversity contri- 
butions and are neglected. Note that the power suppressed admixture of "H in the second odd 
harmonic is consequently neglected, too, and that we also alter the definition of the c^^^^ coeffi- 
cient, since power suppressed twist-three square terms and presumably-small transversity square 
terms are set to zero, cf. formula (2.18) in ^2)- In this approximation the Fourier transform fl^Tl) 
can be evaluated in terms of elementary functions, 

/ I .../,, X \ 
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is considered as an overall normalization factor and the coefficients 
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BH>DVCS C2 
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:^;^DVCS 
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_ (54) 

arise from the higher harmonics of the denominator in fl5T]) . In our case of interest they satisfy 
the condition \a\ > \b\, they are small quantities, and the indicated approximation ( 154|) can 
be considered as justified in a BH dominated regime. Consequently, in such an approximation 
the DVCS-squared term enters only in the overall normalization factor A^. Analogous formulae 
can be also obtained for even harmonics, where we restrict ourselves here to the two lowest 
harmonics. 

Before we discuss the proton case, let us provide a pedagogical example with a spin-zero 
target. In this case we are dealing with three CFFs "H, HcS, and "Ht- Moreover, we assume that 
the second and third harmonics are compatible with zero, which suggest that the CFFs 
and Ht, associated with partonic twist-three and transversity processes, can be neglected. Note 
that this approximation may induce an ambiguity in the phenomenological treatment; instead 
of neglecting "HeS it niay be alternatively justified to neglect the twist-three CFF "Hs, i.e., one 
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sets T-Lcs ~ —xbH- We relate the dominant asymmetries to the twist-two associated CFF "H by 
two hnearized equations 

Ai'^^}'^^ ^ Nc^^n'"^ and A"^""^'^^ ^ Nc^ln''\ (55) 

where we use for later convenience the shorthands H'^"' = 3mn and U^' = DXtH. The coeffi- 
cients 



^^,sin(l</)) 



1 _ 



and Cy,] - 



(56) 

T=0 



T=0 

are calculated from fl22ll24p and fl23|25p and are given as elementary functions of the Fourier 
coefficients s\, c\, and c™. In this procedure, we set the DVCS-squared term in the denominator 
to zero which, however, appears in the normalization factor A^. This overall factor can be 
considered as a bilinear function of the twist-two associated CFF "H or equivalently as a function 
of the asymmetries. To a good approximation, it can be also expressed by the ratio of the BH 
and DVCS cross sections 

r d(j)w((j))daBHi<P) 
j__^d(pw{(p) [daBuW) + daBvcsWl 

The solution of the linearized equations ( 155|) is immediately obtained and, with our twist-two 
dominance assumption, the imaginary and real part of the leading twist-two associated CFF 
reads 



Jmn = ^^Aif,'^^ and ^Rtn = ^^A''°'^'^\ (5^ 
N{A) ^^'^ N{A) ^ ^ 



where Cj^ and Cchc represent now two numbers for a given kinematical point. The normalization 
factor is approximately given as 



1 ..V, ^io,i2 IlJ<Pw{<P)d(^Bycs{<P) 

N ~ r- — — With — 7i — 



jlW 4' ' /:^d0w;(0)rfaBH(0) 



(59) 

'H^a=nT=0 



where /c is a known kinematical factor. Plugging the solution (l58l) into the normalization ([59 
yields a cubic equation in that has two non-trivial solutions: 




iV( A) - - ± ^ 1 - A: {At^T) "-kcl, [a-^'^') ' j . (60) 

The solution with the positive root satisfies the boundary condition N{A = 0) = 1 and is the 
one to take if the BH amplitude overwhelms the DVCS one. Obviously, the solution with the 
negative root satisfies the boundary condition N{A = 0) = and it is the one to take if the 
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unpolarized DVCS cross section is larger than the BH one. Finally, error propagation is done 
in the common manner. If the experimental errors are uncorrelated, we may write the relative 
standard error in the following form 
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where the normalization uncertainty, 
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(63) 



can be easily evaluated by means of (!60l) and one may add a model estimate due to the neglected 
twist-three and other transversity CFFs, which is indicated by the ellipsis. 

The method, i.e., the map of random variables, can be refined by the inclusion of higher 
harmonics and it can be extended to a nucleon target. According to the experimental observables 
and the assumptions, we have generically a set of charge-odd asymmetries, arising from the 
interference of BH and DVCS amplitudes, that we arrange into an m-dimensional vector: 



(A^)^ = (Ai,-- - ,A 



(64) 



For a complete twist-two DVCS off-the-nucleon analysis, we need 4 even and 4 odd harmon- 
ics. Inclusion of the twist-three sector increases the number of harmonics to 16, while for an 
hypothesis-free treatment we need 24 asymmetries. Surely hypotheses can be refined, e.g., in 
[23] the assumption 3m £ = is made in a twist-two analysis. Up to an overall normalization, 
we can formulate a linear mapping problem as follows: 



= N{A^\g) c"^ ■ :f + N{A^\g) b-g 



(65) 



where ^ is the set of sub-CFFs (written as an m-dimensional vector) that one wants to extract 
from the measurements and g is the set of sub-CFFs that is considered as inaccessible within 
a given hypothesis. The inverse m x m coefficient matrix and m x {n — m) matrix h of 
the inhomogeneous term are calculated from the theoretical formulae, e.g., for asymmetries 
proportional to the interference term, we may use 
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Obviously, the solution of the inhomogeneous problem ( 165|) provides us the sub-CFFs in 
dependence on the m observables A"^ and the (n — m) dimensional set of unknown sub-CFFs Q: 

1 



N{A^\g) 



b-g. 



(67) 



The variation of the solution w.r.t. both the measurement and the unknown sub-CFFs can be 
written in terms of an m x m matrix and an m x (n — m) matrix 
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where ® symbolizes the direct product of two vectors, e.g. 
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In the case that (with a given hypothesis) the number of sub-CFFs matches the numbers of 
observables, our equation (1651) reduces to a homogenous one, i.e., we can set the inhomogeneous 
term in fl67|) . containing the matrix 6, to zero and the variation of the solution is given by the 
Jacobian fl68ll. which can be now written in the form 
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(70) 



From this we can easily transform the (variance-) covariance matrix of the experimental mea- 
surements to that of sub-CFFs: 
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where the result reads more explicitly as 



cov(^) 



1 



c ■ cov(A-^) ■ c""" 



(71) 



(72) 



c ■ cov 



{A-) 

dN{A^] 
NdA^ 



dN{A^ 
NdA^ 



^ ■cov(A^) ■ (V® 



dN{A 
NdA^ 



dN{A^] 
NdA^ 



COV 



This representation allows us to discuss separately the overall normalization error given by the 
last term on the r.h.s. of this equation. The normalization factor is determined as in our 
toy example from substituting the resultant sub-CFFs into (l57|l . which provides again a cubic 



21 



equation for A^. The correct solution is picked up by the requirement that is real valued, lies 
in the interval < < 1, and from the knowledge if the observables arise from a BH (A^ > 1/2) 
or DVCS (A^ < 1/2) dominated scenario. If such a solution does not exist, e.g. due to some large 
statistical fluctuation of a mean value, the mapping method is strictly-speaking not applicable; 
however, as we will see below in Sect. I3.3.1l in some of such cases it can be still considered as a 
useful tool. 

If no further experimental information in a given kinematical bin is available, one may use 
a model estimate for the remaining unknown sub-CFFs Q and propagate the estimated uncer- 
tainties by means of a covariance matrix, providing us with the sum 
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(73) 



A more appropriate method would be to constrain the uniformly randomly distributed values of 
unknown sub-CFFs Q and propagate the uncertainties numerically to the errors of the extracted 
sub-CFFs J- ^ formally written as 



dA 




(74) 
(75) 



A third possibility is that the remaining unknown sub-CFFs Q could be extracted from 
asymmetries in the charge odd sector, i.e., related to the DVCS-squared term. Let us arrange 
these asymmetries as an n — m dimensional vector 
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(76) 



and let us suppose that we can complete the number of observables, written now in terms of a 
n dimensional vector 

A^ = {Ai, • • • , A^, A^+i, • • • , An) ■ (77) 

The charge-even observables are linear or quadratic in the remaining sub-CFFs Q, i.e., they are 
constrained by a system of linear constraints. 
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such as those accessed by the transverse target spin asymmetry ( 138|) . or quadratic equations, 

^Dvcs ^ ^ ^T(^) g + g' c^ g], ^ e {m + l, ■ ■ ■ , n} , (79) 
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where the overall normalization can now be considered as a function of the complete set of 
asymmetries (1771) . Substituting the solutions of the linear constraints ( E7^ into ( 175]) and/or ( 179]) . 
one will generally end up with n — m quadratic equations 

Afvcs = jv(^) [a,(A^|iV(A)) + bJ(A^|iV(A)) ■ Q + Q' ■ c,(A^|iV(A)) ■ Q] , (80) 

where all coefficients explicitly depend only on the charge-odd asymmetries and the normal- 
ization factor A^, which appears on the r.h.s. as A^^ with pG{ — 1,0,1}. The quadratic constraints 
fl80p may be analytically solved, too, and the solution may be written for convenience in terms 
of a matrix equation 

a=^c,.cs(A|iV(A)).A--, (81) 

where Cdvcs is an (n — m) x (n — m) matrix that explicitly depends on the asymmetries and on 
the normalization, in general in a non-linear manner (containing roots). Of course, the solution 
(ISTj) is not unique and only real- valued sub-CFFs are of interest. Moreover, one would naturally 
utilize boundary conditions to select the desired solution. Other disadvantages of using charge- 
even rather than charge-odd observables are that the former asymmetries are expected to be 
smaller than the latter ones and the Cdvcs matrix is only known within some given experimental 
uncertainty. The gradient of the remaining sub-CFFs is calculated from 



" iV(A) ''"""^^^^ OA +iV(A) OA ^ ^^N{A)dA' ^^^^ 



where dA°^'^^ /dA is a (n — m) x n matrix that projects on the subspace of charge even asym- 
metries. 

As we have shown, we can map, under certain assumptions, the measured observables into 
the (sub) space of CFFs. We may collect these into the vector 



where the sub-CFFs Q ■ ■ ■ , Qn-m are now a part of J-" and n is now replaced by the number 2n. 
Combining the linear and quadratic solutions (18 ip allows us to find these 2n sub-CFFs from the 
asymmetries by means of 

^=^c(A|Ar(A))-A, (83) 

which can be non-linear in the measured asymmetries A, specified as vectors in ( 1771) . The 
covariance matrix is calculated in the common manner as in (17T|) . replace there A"^ by A, where 
the gradient is now given as a 2n x 2n matrix 

d^F 1 r , dc{A) ^ ^ dN{A)'\ 
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Let us recall that the normalization factor as a function of the 2n asymmetry measurements 
follows from substituting the solution of (155]) into the definition ( IF?]) , which also allows us to 
calculate its gradient: 

dN{A) 



NiA) 
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NdA 

Remaining observables that are not used for the extraction of CFFs can be evaluated from the 
solution ( 183|) and may serve as a test for the validity of assumptions. 



3.3 Local extraction of Compton form factors 

As explained in the preceding section, the extraction of CFFs can be considered as a map 
of random variables from the space of observables to the space of sub-CFFs rather than as a 
fitting problem in which one relies on a given model and tries to find its parameters by means of 
statistical methods. Of course, one may consider the theory as a model and use the least-squares 
method (or maximum likelihood estimation) to extract locally CFFs from measurements. In the 
following sections we use both points of view in local extraction procedures and confront the 
findings. 

The fourteen Hermes measurements, which we will utilize in different variations, together 
with the definition of asymmetries, described in Sect. 13.11 are listed in Tab. [21 Naturally, we will 
arrange the selected asymmetries as a vector, e.g., an eight dimensional one A^ = [Ai, ■ ■ ■ , Ag). 
In our analyses we assume that experimental errors are normally distributed, and the error 
propagation is performed as described in Sect. 13.2. 1[ The correlations of experimental errors 
have not been analyzed to the full extent; however, the results for the overall values as published 
by the Hermes collaboration may suggest that the covariance matrix is mostly diagonal |22j . 
Thus, we may safely assume in the following that the experimental errors are uncorrelated, e.g., 
a eight-dimensional covariance matrix 







cov(A) 












\ 



V 








6^Ay 



^6) 



6'As J 



contains only diagonal entries, which for simplicity (actually, due to missing information we 
cannot do better) are given as the sums of squared statistical and systematic errors. 

As there is no broad consensus for the precise definition of CFFs, we adopt here the con- 
ventions of [16] for unpolarized and longitudinally polarized target asymmetries and take for 
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Observable 


Definition Formulae 


Data from ref. 


Asin{l<j>) 

ACOs{l<j>) 


/",^Alu,i(0) m and 

/" '^''"'''"''Anid) (m and (12711281) 


[22], see Sect. 12:^1 
[22], see Sect. 12.21 


.cos(0(/)) 


/lS^c(0) m and dSEH]) 


[22J, see Sect. |2.2| 


.sm{l(^) 
.cos{l(^) 


/^^^?^AuL,+ (0) (m and P9113TD 

r,^^ALL,+(0) (m and mm 


Tab. 4 from [21] 
Tab. 4 from [21] 


.cos(0(/)) 


nS^LL,+ (0) (m and gTEHD 


Tab. 4 from [21] 


1 sin(ip) cos(l(^) 
^UT,I 

J cos{ip) sin{l(^) 
^UT,I 

A sin{ip) cos{0<j>) 
^UT,DVCS 


J\^^J\^^^A^^j{<P, ^) (USD and (1321341) 
r,^^r,^AuT,i(0, y^) (USD and PWD 
/:.^/:.S^UT,Dvcs(0, y^) (USD and dMl) 


Tab. lb from [T3] 
Tab. lb from [13] 
Tab. la from [13] 


Asin{ip) cos(0<^) 
^UT,I 


/:.^/:.S^uT,i(0, (USD and (nasi 


Tab. lb from [T3] 


Asin{ip) sin(l(^) 

1 cos(ip) cos(l<7i) 
^LT,I 

A cos{ip) cos(0(/>) 
^LT,BH+DVCS 


J\^^r/-^A^^j{^, ^) m and (1331341) 
/",^^^r,^^ALT.i(0, V^) (HID and (1361137D 

/:.^/lS^LT,even(0, (ffl and dSHD 


Tab. 2 from [T5] 
Tab. 2 from [15] 
Tab. 3 from [I5] 


^ cos{ip) cos{0<f)) 
^LT,I 


/:.^/:.S^LT,i(0, ^) dlZD and (IISETD 


Tab. 2 from [15] 



Table 2: Observables from Hermes measurements that are utilized for the extraction of twist- two 
associated CFFs. 



transversally polarized ones the 1/Q expanded expressions from [14]. Moreover, we will restrict 
ourselves here to the twist-two sector, i.e. to the CFFs J-" G {"H, "H, related to the observ- 
ables by setting the remaining eight CFFs J^cs and J-t to zero, which is justified by the fact that 
higher harmonics are compatible with zero or are difficult to interpreto. The imaginary and real 
parts of the twist-two associated CFFs J-" are collected in an eight dimensional vector 



7im cdm 'c-"^ 



cJm c njint njinc oytc c 



(87) 



^In contrast to other higher harmonics, in the longitudinal single target asymmetry a relatively large sin(20) 
moment with large uncertainties has been observed, which increases with growing ~t and which is maximally 
up to ^ 2(7 deviations away from zero |21j . However, this large value is concentrated in a single bin in the xb 
projection, which is supportive of concluding that the amplitude is simply a fluctuation. 



25 



written below as T 



in terms of two four dimensional column vectors 



J= 



■3m 



(n\ 
u 

£ 



and ^ 



•IHc 



£ 



where we use £ 



Xb 



-£. 



The (approximate) redefinition of £ into £ removes the factor xb/(2 — xb) that enters in the 
common form factor decomposition and makes the kinematical coefficients of £ in the cross 
section rather small. This redefinition restores also the common "Regge" behavior, i.e., £ oc 
Xg"^*"* for small xb (see [13]) and simplifies the discussion of the real photon limit for CFFs [T7] . 
However, the common definition is better suited to relate GPDs to form factors. 



3.3.1 Maps of asymmetries to CFFs 

To determine the sub-CFFs flHHjl by means of a map fl83|l . we select eight twist-two related asym- 
metries out of the measured first harmonics, listed in Tab. |21 Thereby, we take the kinematical 
means from Tab. [U In any of the possible maps we assume the twist-two dominance hypotheses 
and we naturally employ for the access of the imaginary parts of the twist-two associated CFFs 
J-" (specified in (l88l) ). the dominant single spin asymmetries in the following sequence 

n 

£ 



A' 







/ .sin(10) \ 
^LU,I 


A2 




.sin(10) 


^3 




^ sin(i/j) cos(10) 
^UT,I 


^4^^ 




/,cos(</3)sin(10) 
\^UT,I / 



-3m 



3m 



(89) 



see the expressions fl26|29|32|[35|) for asymmetries in terms of CFF combinations ( l28|31|34|37p . 
The analogous relation for the corresponding even harmonics and the real parts of CFFs reads 



A-r 



( 



.C0S(1(/)) 

.cos{l(^) 

sin(</3) sin(l(/!)) 
LT,I 

< COs(</3) COs(l</>) 

V^LT.I 



A"- 



-We 



u 

£ 



(90) 



see f l27|30|33|l36|) and f l28|31|34|[37|) . In the following we utilize, however, two different analytic 
methods for the extraction of the real parts. 

First, we consider a linearized map ( 1671) . where we take only charge odd asymmetries that 
arise from the interference of BH and DVCS processes. Thus, we transform the single longitudi- 
nally polarized target spin asymmetry J^^^f^ and the longitudinally double spin-flip asymmetry 
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+''^'* (both measured only with a positron beam), appearing in the four dimensional vectors 
fl89|) and fl90|) . into the charge odd sector. In doing so, we eliminate the interference term in the 
denominator of these asymmetries by means of the known beam charge asymmetry, i.e. we find 
from the relation ffTOl): 



<sm(10) 
"^UL,I 


.sin(10) 

~ 4sin(l9i) 
~ ^UL,+ 


1 _|_ y^cos(0(/)) 

X _|_ y^«5s(0(^) 


1 ^cos(24,) 

2 ^ 


,cos{l<f>) 


acos{14,) 


1 _|_ ^cos{0(/)) 


l^cos(20) 
2 ^ 




ACOs{l<f>) .C0S{1(^) 

"1 ■ ^LL,BH ^LL.DVCS ' 








.cos(l</>) 
^LL,BH 



sin{2(^) 
UL,+ 



a: 



cos(l(/!)) 
C 



A 



cos (30) 
C 



+ 



sin(10) 
UL,DVCS ' 



(91) 



, <cos(0(/!>) ^cos(l(/)) 



2 



cos{2(^) 



^cos(l(/>) ^ ^cos(3(/!>) 



(92) 



cos(2(/!>) <sin(20) tcos{2<j>) ^sin(l(/!>) 



4 cos(l(/!>) 



^UL,DVCS' ^LL,DVCS' 



and we neglect the twist-three related asymmetries Afj , , 
and higher harmonics, indicated by ellipses. Moreover, for the charge odd double longitudinal 
spin asymmetry, the first harmonic of the BH contribution turns out to be small and can thus 
be safely subtracted (proton form factor uncertainties are neglected). The imaginary and real 
parts of the CFFs follow from the solution of two linearized homogenous equation, see f E7|) with 
6 = 0, which we write in the form of the map fl55]) as 



N{A) \ 



4x4 




(93) 



Here 0„xm denotes a n-by-m zero matrix in the 8x8 coefficient matrix c, calculated by the 
inversion of fl5^ . and the first odd and even harmonics 







/ .sin(l<^) \ 
^LU,I 


A2 




.sin(l(/)) 
^UL,I 


As 




< sin(i/5) cos(l(/i) 
^UT,I 


[aJ 




.cos((/3)sin{l(/)) 
\"^UT,I / 



and A"" 



fA,\ 

A, 
A, 



/ A 



cos(l(/)) 
C 

cos(l(^) 



\ 



4 

^LL,I 
.sin(i^)sm{l(/)) 
^LT,I 

< COs(ly3) COs(l(/)) 

\"^LT,I 



(94) 



respectively. Consequently, in this linearized map, the cross talk of imaginary and real parts of 
sub-CFFs arises only via the overall normalization factor, which we analytically determine from 
the cubic equation (!57|) . The gradient of this linearized sub-CFFs solution can be analytically 
calculated by means of (!70|) and the covariance matrix follows from ( !7T|72|) and ( 186|) . 

Unfortunately, the statistics in the double spin-fiip asymmetry measurements are rather 



limited, in particular, for the ^lt i 



sin((/3) sin{l(^) 



and A 



cos(i^) cos(l(/!>) 
LT,I 



harmonics that essentially constrain 



27 



the real value of £ and £. Hence, we try a more general map (l83i) where we replace the two 
longitudinal-transverse double-flip asymmetries in the charge odd sector by the two DVCS- 
squared related asymmetries ^ut dvcs'^'^'' ^^"^ ^ll +''^^ which are in general measured with smaller 
uncertainties, exemplified in Fig. [2l We arrange our eight asymmetries now into three parts. 



With = f = f ^ ... , ^ = '^UT.DVCS 




.cos(l(^) / ' \ A ] \ ^cos(0</)) 



(95) 



and is the same vector as in flMjl . For the first even harmonics of the charge odd asymmetries 
A^°^ we take here the first two entries of the corresponding four dimensional vector in f lM|) . 
which are mostly sensitive to the real parts of CFFs "H and "H. Consequently, we split the 
four-dimensional vector J-^^ into two two-dimensional vectors. 



= 1^ ^"^^ j with n^' = ^He j and = ^He | ~ ) . (96) 

Here, "H^' (= ^ in the notation of Sect. 13.2. ip is considered as a solution of the inhomogeneous 
equation ( 167|) and £^'^ (= Q in the notation of Sect. 13.2. ip governs its inhomogeneous term, i.e., 

n""' = ■ A=°^ - c^, ■ b^, ■ £:«^ . (97) 

The transverse single target spin asymmetry in the charge even sector is proportional to the real 
parts of CFFs £ and £, see (l38l) . and provides us one linear constraint 



Aj = A'^!;,%™f = N{A) cl{A'"\ A^°^) ■ £:^^ (98) 

where the cj-coefficient depends on the six charge odd asymmetries A^^^ and A'^°^ or sub-CFFs 
_^3m 'j-lp^^ that are extracted from the set of linear equations. Our system of equations is 
completed by the measurements of the longitudinal double spin flip asymmetry (H7p . 

A <cos(O0) <cos(0</<) I <cos(0</)) I tcos{0<j)) 

^8 = ^LL,+ ~ ^LL,BH + ^LL,I + ^LL,DVCS • 

It contains (in addition to a large BH-squared term) contributions from the polarized interference 
term and the twist-two DVCS-squared contribution in the numerator. Here we cannot assume 
that the interference term overwhelms the DVCS-squared term; rather we assume that the two 
terms enter kinematically on the same level. Hence we have a quadratic constraint for the CFF, 
which we write as 

As = = N{A) [a(A^'^, A^°') + W{A'"\ A^°^) ■ £^' + £^'' ■ c{A'"\ A^°') ■ S'^'] , (99) 
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where all coefficients depend again on the six measurements A*''" and A^"^, used previously. The 
linear constraint together with the quadratical one fl99p can be analytically solved, giving 
us the desired last two sub-CFFs which we write as a matrix equation, too 

1 



N{A) 



ce{A\NiA)) ■ A^ 



(100) 



where the 2x2 matrix C£{A\N{A)) is a non-linear function of the eight asymmetries A. The 
final solution can be written as in ( 183|) with a 8 x 8 coefficient matrix c{A\N{A)), 



N{A) 



04x2 
02x4 C(Rc 
02x4 02x2 



-c^, ■ b^Hc ■ ce{A\N{A)) 
ceiA\NiA)) 



I 



\ 



^cos 
^DVCS 



(101) 



The normalization N{A) is again determined by the consistency equation f l57l) . The gradient 
of the sub-CFFs solution is calculated by means of (!84|) and the covariance matrix cov(jF) is 
calculated from (ITTIl . replaced there A^ by A, and (186|1 . 
A few comments are in order. 

It is clear that, within our hypothesis of twist-two dominance, the matrix valued equations 
fl93|101l) are solutions of well defined problems. However, in practice we have to deal with the 
fact that the eigenvalues of (or its inverse c) matrix can become very small (big). Obviously, 
this unpleasant property arises from the fact that in particular the helicity target flip CFFs £ 
and £ are kinematically suppressed in the set of observables. 

In addition we have six observables, namely the three lowest harmonics fl41f42f43p 



A 



cos{0<j>) 



A. 



sin{i^) cos(0(/!>) 
UT,I 



A 



cos(ip) cos{0(/)) 
LT,I 



and three twist-two related asymmetries fl38|39|47p 



AS\n(tp) cos(Ocf>) < cos((/3) cos(0(/)) ^cos(0(/)) 

^UT.DVCS ' ^LT,BH+DVCS ' ^LL,+ 



for our linearized map ( 193|) . 



or alternatively the asymmetries fl33|36|39p 



A 



sin(i^) sin(l(/!>) 
LT,I 



A 



cos{ip) cos(l(/i) 
LT,I 



^LT?BH+Dvcs fo'^ our map (JM]) 



which are predicted by the extracted twist-two associated CFFs and serve as a consistency check 
of our extraction procedure. However, we should bear in mind that the three lowest harmonics 
fl41|42|43p may be more strongly contaminated by twist-three related CFFs. 

For comparison we will also employ the brute force method, where we solve eight quadratic 
equations numerically and numerically evaluate the variation of the solution. This allows us 
to judge the validity of the approximation used to linearize the constraints. Thereby, we can 
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employ the original Hermes data for the first even and odd harmonics ( I891I90I) rather than 
the transformed ones f l94p . As for the linearized map, we obviously have six observables for 
consistency checks available. 

Let us now present our results. Within the linearized map we find that the consistency 
equation for the normalization fl57j) has two non-trivial real-valued solutions in eleven out of 
the twelve kinematic bins and the overall asymmetry values, shown in Fig. O According to the 
experimental indications, we pick the solution for the BH regime, where we find for the weighted 
cross section ratio 



0.75 < 



J^^d(l)w{(f))daBH{<P) 



< 0.95 



(102) 



The inverse transformation of the sub-CFF solution back to the original asymmetries, where we 
use the set of non-linear equations, reproduces in general the means and standard errors of the 
seven asymmetries 

J <sin(l(/>) <sin(l(/)) < sin(iy3) cos(l<)!)) »cos{ip) sin(l<j>) ACOs{l<j>) Acos{0(f>) Asin{ifi) cos{0<f>)\ 
|^LU,I ' ^UL,+ ' ^UT,I ' ^UT,I ' ' ^LL,+ ' ^UT,DVCS J ' 

(the derivatives of which we have used for the map to the sub-CFFs space) in most of the cases 
on the level of a few percent. We observe in some bins a larger deviation only for A^-[^^^\ 
which is naturally explained by the fact that we neglected in our linearization procedure the 
DVCS-squared term, which becomes important if the asymmetry is small. However, also in this 
case the original data can be considered as well-reproduced. This inverse map is shown for all 
considered observables in Figs. [814TT1 shown in Appendix El as empty circles. We add that a 
typical mapping example is presented for the overall asymmetry values in Fig. HJ 

A true one-to-one map of random numbers is reached if we use the brute force method, 
where the start values may be taken from the solution of the linearized map. Moreover, the 
mean values of the remaining asymmetries, used for the consistency check, are well reproduced 
on the ~ la level, see stars in Figs. [3] and Hi 

Let us have a closer look to the problematic bin ^3. The linear map (p3|) yields 



Jm 



n 
s 

V ^ J 



N{A) 



( 7.6 ± 1.2 \ 
1.8 ± 1.3 
-4.5 ±6.0 
\ 11.0 ±6.1 / 



and SHe 



( U\ 

n 
s 

\ ^ J 



N{A) 



( 0.4 ±2.1 \ 
5.7±4.1 
-28.7 ±21.2 
\ -24.8 ±22.2 / 



,(103) 



where the uncertainties are calculated from ( 1721) by neglecting the variation of the normalization. 
The consistency equation (1571) provides us two complex valued solutions = 0.5 ± 0.06 i, where 
the smallness of the imaginary part tells us that the inconsistency is rather weak. From (11031) we 
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may conclude that the large mean values for £^'^ and , which suffer from a large uncertainty, 
ruin in turn the normalization constraint (1571) . Indeed, setting one (both) of them to zero allows 
us to solve the normalization constraint, where ~ 0.75(0.85). There are other possibilities 
to turn around the normalization inconsistency. In our case the data mapping and/or the 
linearization procedure imply this inconsistency. Applying the brute-force method to the eight 
original observables provides the solution: 



£ 



11.8 ±8.9 \ 








2.6 ±2.4 










and 




-8.1 ± 13.2 




£ 




15.6 ±12.7 ) 









( 



0.8 ±3.6 \ 
6.2 ±9.9 
-41.9 ±52.4 



(104) 



where the errors are propagated as discussed above. From these findings we obtain the relatively 
small normalization factoio = 0.64 ± 0.28. Substituting this number into (I103P shows us that 
the sub-CFF mean values from the linear map are compatible with those of (I104p . obtained with 
the brute-force method. However, in particular, the net error for the sub-CFF Jm'H, calculated 
from the sub-CFF errors (11031) and the normalization uncertainty, turns out to be smaller than 
that in (fTOij) . 

The resulting CFFs from the linearized (circles) and one-to-one (stars) maps are presented 
for all twelve bins in Fig. El Both maps provide rather similar results, except that the errors 
for the sub-CFF Jm'H in bins ^^3, #7, and ^11 are approximately two times larger. In these 
three bins the cross section ratio < 0.7, calculated from the one-the-one map, is rather small. 
We conclude that non-linear effects in these circumstances are rather important for the error 
propagation (see above discussion for bin #3). Clearly, the imaginary part of 1-i is (as expected) 
positive, rather large and incompatible with zero. Note that the relatively large errors in bin 
#5 and #9 are (partially) induced by the fluctuations of 3x^£] see for example the analytic 
expressions for the beam spin fl26ll28p and the cos0 projection ( 132|M1) of the single transverse 
target spin asymmetries. 

All other sub-CFFs can be considered within the uncertainties as compatible with zero. Sur- 
prisingly, the imaginary part of 1-i possesses even smaller absolute errors than ^mH. Both of 
these sub-CFFs are contaminated in the first place by proton helicity flip sub-CFFs 3m £ and 
Jm£^, respectively (see CFF combinations fl2S]) and (EI])), which both suffer from larger uncer- 
tainties. However, comparing the approximated expressions fl2B]) and fl2^ for the asymmetries 



■^We neglected here correlation of sub-CFF errors otherwise the normalization error increases to 0.48. In turn 
we also neglect the error correlation in utilizing this normalization error in the linear map ()103|) . 
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/|Sin(l</') „„J /isin(l</.) 
^LU,I ^""^ ^UL,+ ' 

AtvT 1/(2-1/) Jm[Fi?^-^F,g + ---] 

one realizes that the beam spin asymmetry has an additional relative suppression factor, 

~ 0.5 with xb ~ 0.1 and ~ 2.2 GeV^ 



for typical Hermes kinematics. Hence, the larger errors of the polarized longitudinal target spin 
asymmetry are reduced in the propagation to the sub-CFFs by a factor of two or so. Another 
reason why the absolute error of JmH is larger than of JmH is that the former sub-CFF is 
more sizeable and so it is also more important for the normalization ratio f l57p than the latter. 
In return, JmH suffers from a larger absolute error, see, e.g. the explicit form of the covariance 
matrix ( !72|) for a linearized map. Note, however, that the twist-two hypothesis may induce here 
a bias, since we neglected here the twist-three induced siiKp harmonic of the DVCS-squared 
term. 

The imaginary parts of the proton helicity flip CFFs S and S are less constrained, partially 
due to larger errors of the single transversely polarized target spin-flip asymmetries — however, 
also in part due to their being kinematically suppressed, see ( 132| |34|) and fl35|37p . Only the real 
part of the CFF "H, also compatible with zero, can be definitely considered as well constrained 
due to the precise data and since "H dominantly enters in the charge asymmetry fl27|28p even 
without additional y suppression. As one sees, the real parts of the remaining three sub-CFFs 
are very noisy, in particular, those of the proton helicity flip CFFs. A generic GPD model 
interpretation of our sub-CFF findings is presented below in Sect. 13.41 

In the alternative map we replace the charge odd asymmetries ^^^^^'f'^ and y^^^(|^)'=°'^(i<^) 

by j^^^^^'l'^ and ^ut dvcs°'^'' ^^^^ contain the DVCS-squared term in the numerator. The nor- 
malization is now determined by the non-linear equation, see also ( 159|) . 

N = N{A, N) with N eR and < < 1 , 

the solution of which yields now four roots. Literally taken, the mapping method is now only 
applicable in 7 out of 12 bins, namely, in 

#1, #2, #4, #9, #10, #11, #12, 

and the overall bin. The failure of the method in bins #3 and #6 is caused by the longitudinal 
double spin asymmetry measurements A'^-^^^'^\ which as mentioned above contradicts the as- 
sumptions that the BH amplitude overwhelms the DVCS one, see Fig. [31 The failure in bin #8 
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Figure 5: Resulting sub-CFFs from a linearized (circles, shifted to the left) and a one-to-one map (stars) 
of eight twist-two dominated charge odd asymmetries as well as from a least squares fit (triangles, shifted 
to the right) to fourteen twist-two related observables for each of 12 Hermes bins. 

is related to the large transverse spin asymmetry ^ut dvcs"''^'' (^^^ -^^S- SI) while the inconsisten- 
cies in the mapping method for the bins #5 and #6 are rather weak. Hence, we conclude that, 
even if the uncertainties of the both replaced charge odd asymmetries are rather large, these 
observables yield important constraints that ensure the consistency of the mapping method. 



3.3.2 Local least squares fits to asymmetries 

Let us now employ the method of least squares, where one looks for the minima of the function 



A - a{:f) 



A - AiJ^) 



(105) 



Here, the vector A contains the measured asymmetries and cov ^(A) is the inverse of the 
covariance matrix. For uncorrelated errors the covariance matrix ( ISBjl and its inverse are diagonal 
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and the x^-function (11051) reduces to the most common form 



X 



A, - AAJ=\ 



which is well-known. The error propagation is often performed via the Hessian matrix 

/ i^ii ■ ■ ■ Hxn \ 



H 



1 



with Hi. 



J 



2 dJ^idJ^j 



(106) 



X Xn 



y Hjil ' ' ' Hjin 

where its inverse provides the covariance matrix for the sub-CFFs 

covCF) = H K 



(107) 



Let us first remind the reader that, instead of finding a one-to-one map of normally distributed 
random variables with the methods outlined in Sec. 13.2.11 and used in Sect. 13.3. Ij one may 
equivalently utilize the least squares method. Obviously, if a solution (I83p exists, written as 

:f = :f(a) with a = a(:f), 

the function (llOSp takes exactly the value = 0- However, using blindly a "black box" 
fitting routine to extract the CFF values from DVCS asymmetries, one may find only one global 
minimum with ~ due to numerical errors. This solution may be associated to an unphysical 
root, e.g., in Hermes kinematics it may be associated with a solution where the DVCS cross- 
section dominates the BH cross-section. Hence, in such "fits" one must search for all ^ ^ 
minima and decide then by means of the cross section ratio (1571) which of them is to be considered 
as the physical one. Alternatively, one may implement the constraint N{j^) > 1/2 for the BH 
dominated regime N{jP) < 1/2 for the DVCS dominated regime or one can constrain the 
value of a single sub-CFF that is not very well determined by the data. It is easy to realize from 
(I105H107]) that, if a one-to-one map exists, the covariance matrix (I107p can be written in the 
form 



cov(^) 



9:f 



-1 



cov(A) 



-\ -1 



(108) 



or by means of the inverse function theorem in the form of fl7T|) . constructed from the Jacobian 
flH^ . We emphasize that the use of constraints (particularly if the fitted parameter ends on the 
boundary) may influence the results and so a true one-to-one map cannot be obtained. 

We used the equivalence of brute-force and least squares methods, employed for the set of 
charge odd asymmetries (IMIl . for a numerical cross check between two independent software 
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tools. Utilizing the popular minimization routine MINUIT [H] in one code and the brute-force 
method in the other, we obtain (except for bins #3 and #4) the same one-to-one map that is 
shown in Fig. |5] (stars). In these two bins minuit finds global minima in the DVCS dominated 
regime rather than the BH one, exemplifying that taking here the global minima gives an answer 
that we consider to be wrong. The small deviation of the total 

= 2.95 X 10^^ or xV^d.o.f. = xV(12 x 8 - 8) = 3.35 x 10"^ 

value for all 12 bins from zero is here to be considered as a measure of the numerical accuracy 
rather than reflection of statistical fluctuations. 

Obviously, we can employ the least squares method to an ill-posed mapping problem, e.g. as 
discussed in Sec. 13.3.11 for the alternative map with the set of eight observables (p5l) . where an 
inconsistency appeared in a few bins. Of course, in such a case the value will differ from zero 
and can be taken as a measure of the inconsistency, which may originate from the statistical 
fluctuation of means in data. 

Moreover, we can use the least squares method also in the case of an overcomplete set of 
equations. In such a case the value of a global minima will also differ from zero. Including the 
two twist-two associated observables ^ut dvcs''''^^ ^lt ra+ovcs interference 
term dominated asymmetries fl89|90l) . the total = 14.3 value differs now significantly from 
zero. Note that 

xV^d.o.f. = 14.3/(12 X 10 - 8) ^ 0.13 

is naturally a small number and that the value in a given bin can be considered as a measure 
of the size of statistical fluctuations and/or the validity of the utilized twist-two dominance 
hypothesis. The resulting sub-CFFs are entirely compatible with those from the one-to-one map 
(not shown). 

Extending the twist-two dominance hypotheses also to the lowest asymmetry harmonics, we 
have 14 observables available and we may consider the twist-two associated sub-CFFs as eight 
separate independent parameters that we would like to extract. Formulating the extraction 
problem in such a manner makes closer contact to the work in [25], where all 24 asymmetries 
that were available at that time were used. In our fit, the total value increases to 

X^ = 67.7 or xV^-d.o.f. = 67.7/(12 x 14 - 8) ^ 0.42 , 

the mean values remain in general stable (see triangles in Fig. [5]) except for bin 7^8 and #12, 
in which the solution in the DVCS dominated regime is again obtained. For 7^8 we present the 
MINUIT outcome and one realizes that the solution in the DVCS dominated regime possesses 
a very large uncertainty. For bin 7^12 we used the constraint I^Hei"! < 40, to find the local 
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minimum that provides the BH dominated solution. Thereby, the fit ends on a boundary 9^e E = 
—40, and, hence, the error propagation for yicS completely fails — see the corresponding triangle 
in bin 7^12 on Fig. [5l Compared to the one-to-one map, the uncertainties naturally decrease 
and become even smaller than in our linearized map (filled circles). As a consequence of this 
decrease, one may view now the real part of S as not any more entirely compatible with zero. In 
conclusion, extending the twist-two dominance hypothesis also to observables that are potentially 
more contaminated by the remaining sixteen sub-CFFs and/or the inclusion of constraints may 
underestimate the errors and can potentially result in an overinterpretation of the resulting 
sub-CFFs. 

Finally, let us compare our findings with those obtained from a regression analysis. Here, 
one would consider the value as a statistical measure and would consider a solution with 
the value of x^/d.o.f. ^ 1 as optimal. As we have seen, the precision of data and the fact that 
for any given observable most of the sub-CFFs are kinematically suppressed, presently prevents 
us from accurately determining all eight sub-CFFs. Although we have already answered the 
question 'Which sub-CFFs can be extracted with some reliability?', in the rest of this section we 
attempt to address this same question using the method of stepwise regression. Thereby, one first 
performs eight separate single-sub- CFF fits and sees which sub-CFF alone describes the data 
best (measured by value). Then one proceeds to two-sub-CFF fits, where the best sub-CFF 
from the first step is combined with each of the seven remaining sub-CFFs and the best-fitting 
pair of sub-CFFs is retained. This procedure is continued until there is either no improvement in 
the description of the data or new sub-CFFs are not extracted with any statistical significance. 

To make this stepwise procedure more reliable we temporarily remove from consideration the 
observables ^ut dvcs^'^'' ^^"^ ^ur m+dvcs because, for these harmonics, squared combinations 
of sub-CFFs play a dominant role so a) they cannot be reasonably described in first-step single- 
sub-CFF fits and b) their inclusion would introduce strong correlations between sub-CFFs in 
the second and further steps thus potentially introducing bias for some sub-CFFs. 

Fitting was performed by standard minimization of the function separately in each kine- 
matic bin, using the minuit package. As expected, it is JmH that gives clearly the best 
description of all data. The single-sub-CFF fit of JmH gives /''^d.os. = 198.4/132, with second 
best being mzTi with xV^d.o.f. = 472.7/132. 

In the second step, we have two equally good fits: 

• fit of Jxnn and ^cU with xV^d.o.f. = 102.3/120, and 

• fit of 3mn and 9^c£ with xV^d.o.f. = 103.0/120. 

(Next best being fits to 3mn and ^HeH with xV^d.o.f. = 122.4/120 and to JmH and ^Hc^with 
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Figure 6: Results of least-squares fits in two scenarios with only a small number of CFFs locally fitted 
to data separately for each of 12 Hermes bins. First, with only 3mT-L and 9\t£ fitted (red diamonds) 
and, second, with JmTi, d\z'H and ^mH (purple pluses). For comparison, result of one-to-one mapping 
procedure from section 13.3.11 is also shown (green stars) . 

}d/^A.o.i. = 185.4/120.) Trying now the third step, for several choices of the third sub-CFF 
quahty of the fits improve somewhat (measured by their p-value) but the values of this third 
sub-CFF cannot be extracted with any statistical significance. Furthermore, after adding the 
^UT.DVcs'''^'' data, the fits of JmH and DlcS improve much more than the other scenarios due 
to the dominant contribution of terms involving the JmH ■ D\tS product to this observable, 
cf. (l38l) . Thus, as our final results we present two scenarios: one selected by this stepwise 
regression procedure and another, more in agreement with common expectations, where fits are 
done with the predetermined set of sub-CFFs JmH, D^e'H and Jm'H. In both scenarios we now 
perform fits to the complete set of 14 Hermes observables. 

• Scenario 1: Fit of JmH and ^icS. x^/''^d.o.f. = 134.2/144, when adding values for all 
12 bins. Fits are bad for bin #3 (xV^d.o.f. = 19.9/12) and bin #8 (xV^d.o.f. = 21.1/12). 
Other bins are fine. 
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• Scenario 2: Fit of 3mH, D^cH and 3mH. xV^d.o.f. = 148.8/144. Bins #3 and #8 again 
show a bad fit result, with x^/^d.o.i. = 21.8/11 and x^/rid.o.f. = 21.1/11. 

We note that, in both scenarios, the fit in bin #8 is fine if the ^ut dvcs''''^'* ^lt ra+Dvcs 
data points are removed. The resulting sub-CFFs from both scenarios are plotted on Fig. |6] and 
are compared with the one-to-one map. One realizes that the sub-CFF 3xnT-L is quite robust 
and that the real part of OFF S differs on some 2a level from zero. The results from the second 
scenario essentially agree with those in [23], showed there for three selected bins #2, ^^3, and 
7^4. We add that those results were obtained by means of model-independent least squares fits 
to twenty-three asymmetries. However, most of the asymmetries were not related to twist-two 
dominated quantities and the sub-CFF parameters were deliberately chosen with reference to 
GPD model constraints. How successfully fits from our two scenarios (and a global world fit 
presented in Sect. 13. 4p describe particular observables is visible on Figs. IHHTTlin Appendix lAl 

3.4 Uses of HERMES data for model builders and in global fits 

Let us discuss the constraints from Hermes DVCS data, presented in terms of CFFs on Fig. |5l 
for GPD model builders. We consider here the GPD framework in the perturbative leading order 
(LO) approximation and we can safely restrict ourselves to the qualitative aspects. We adopt 
and refine here some older discussions, given in Sect. 5.1 of [15] and illuminated with BMK 
model predictions in Sect. 5.2 there. The imaginary part of CFFs is, in this approximation, 
given by the GPDs on the cross over line: 



F{x,x,t,fi^ = Q2) ^i* lamJ^(xB,t, 



As known, the valence quark part of GPD H is essentially governed by the asymmetry A^y^'^^ 
which is almost saturated by the valence quark content of the forward PDF, decorated with 
some t-dependence. On the other hand, rather generic model estimates tell us that sea quark 
contributions are important in Hermes kinematics, i.e. the partonic decompositioij^ reads 

3mn = 3m n^^i + 3m ?/sca with 3m ?/sca ~ 3m n^^i . (109) 

Consequently, the Hermes data, taken in terms of Jm'H, require a small skewness effect for GPD 
H, which is also required for the LO description of Hera collider data. To get an easy handle 



^Our decompositions of CFFs, which are charge even, in valence and sea quark parts contain squared quark 
charges and the terms valence and sea are adopted from the common terminology as used in global parton 
distribution function fits. 
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on the real part of CFF one may use instead of the LO convolution formula a signature-even 
GPD dispersion relation, see [321 ESI EH SS] and references therein 

9^e|^l(a;B,t,Q') = PV rrfx-^|^l(x,x,t,/x2 = Q2)±I?(t,Q2) . (no) 



For Hermes kinematics, we can take for granted that the sign of the resulting real part is 
determined by the "Regge" -behaviour, inherited from the PDF behavior. Hence, valence quarks 
provide a large positive real part while sea quarks contribute a negative part and, in addition, 
there is a subtraction constant that is related to the "D-term" |17]. The Hermes data, taken 
in the form of IHe H, tell us that the modulus of this quantity is rather small and, hence, we 
would interpret it as representing a cancelation between the three contributions. Note that the 
experimental constraint on the subtraction constant depends also on details of the GPD model 
and it is therefore rather weak. 

Going along the same line, we can now discuss the CFF "H. Its real part arises from a 
signature-odd GPD dispersion relation. 




where no subtraction is needed. From phenomenological PDF parameterizations, one expects 
that the GPD H, taken on the cross-over line and as it enters in the DVCS amplitude, is (much) 
smaller than GPD H. Hence, one expects that both the real and imaginary parts of CFF "H are 
relatively small. This is entirely compatible with Hermes data, which surprisingly provide us 
also with a rather strong constraint for JraT-t. 

A generic discussion can be also given for the CFF S. However, here only form factor 
information is available. Hence, one would assume that the zero-skewness GPD has a simple 
functional form in which nodes are absent that arises from Regge and large-x arguments. The 
normalization of such valence quark GPDs is adopted from the anomalous magnetic moments 

K^^^j = 1.673 and k^^^j = —2.033 , 

taken from the nucleon. A crude approximation of the size of the valence part of Jm S is obtained 
if we assume that the same functional form holds for all valence quark contributions, 

3^ ' ^ Jm£:.,~-Jm7^, 

where we used (11091) . However, we should keep in mind that the functional form w.r.t. both 
the t- and x-dependences may alter our estimates. It is expected that the t- and x-fall off for 
E GPDs is steeper than for H GPDs |18], i.e. GPD-model-refined estimates would give even 
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smaller predictions for JraS^g^. Considering the CFF data in Fig. [5l one realizes that the noise 
of 3m S is of the order of the means of Jm'H. Hence, contrary to GPD-model-based claims, e.g. 
that Hermes data provide a constrain on the quark orbital angular momentum decomposition 
[T3] or the suggestion that negative sea quark contributions to 3m. S are favored |19] , our generic 
arguments tell us that a partonic interpretation of Hermes data in this specific case is entirely 
biased by model assumptions. Let us add that the dispersion relation f lllOp that we used for an 
estimate of DltS, together with the standard model assumptions, tells us that the large negative 
^tS scenario that we obtained from the regression method is difficult to understand from the 
GPD model perspective (positive subtraction constant for negative D-teim, positive contribution 
for valence quarks, and positive contribution for negative sea quarks). 

To complete our short examination of the CFF data, we mention that the CFF S should 
contain a pion pole contribution that should be large at small —t [501 El]- Sometimes, GPD 
model builders believe that this is the most important contribution and neglect for that reason 
the imaginary part. If one wishes, one can see sizeable and negative S sub-CFFs at small —t — 
see the first bin in Fig. O where the significance for the imaginary part is even more pronounced 
than for the real part. This clearly contradicts the common GPD model assertion. However, 
since these data are very noisy, a definite conclusion cannot be drawn and we consider the CFF 
S as essentially unconstrained from Hermes DVCS data. 

For comparison purposes, we also performed one global model fit to the world DVCS off-the- 
proton data with a version of the hybrid model used in [38] that we initially used to access the 
GPD H from unpolarized proton DVCS data. The hybrid model comprises a full GPD model 
in the fiavor singlet sector (dominated at small x-q by sea quark and gluon contributions), while 
in the fiavor non-singlet (or valence quark) sector dispersion relations are used and, hence, only 
the GPDs on the cross-over line are needed. Keeping in mind that, apart from Hermes, no 
other experiment with a proton target could provide information on the full separation of the 
various CFF contributions, we neglect in our fit the 3m S and 3m £ contributions, which are 
compatible with zero for Hermes kinematics. However we do include DltS and D\c£, which 
are related to subtraction constants in the dispersion relation f lllOp and in the oversubtracted 
analog of fillip , respectively. The reason for doing so is that HALL A unpolarized cross section 
measurements [52] indicate a rather large unpolarized DVCS cross section, which suggest that 
the DVCS amplitude contains a large real part. In our hybrid model, the contributions of sea 
quarks and gluons are modelled using conformal moments of GPDs H^*^^ and H'^, respectively, 
and LO QCD evolution is taken into account, while the contributions of valence quarks are 
described by directly modelling 3m'H and 3m'H and using dispersion relations to obtain D^e'H 
and fRtT-L (and evolution is neglected). Details of the model are given in ref. [3B], to which we 
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refer the interested reader. 

In particular, we used the following data sets: 

. The Hermes combined data on Ai^^}'^\ A™^^'''^) and A'"""^'^^ [22]; on aJl^'/^ and A™',^^ 
[21] : and on work with statistically independent data we considered 

only the projection of the data along the —t axis, i.e. just the first third of the published 
18 [22] or 12 ^ EI] kinematic bins. This gives 3x6 + 2x4 + 4 = 30 data points. 

• The first (sin</)) harmonics of the CLAS data on a) the beam spin asymmetry with an 
unpolarized target [53], where we used only data with > 2 GeV (4 points), and b) the 
longitudinal target spin asymmetry with a polarized target [51] (6 points). 

• Fourier transforms of the Hall A measurements of beam spin difference (12 points) and 
beam spin sum (8 points) [52] , where cross-sections were weighted with the inverse product 
of the Bethe-Heitler propagators. 

• Measurements by the HI collaboration of 

— The DVCS cross section differential in t, [55], Table 1, 1996-1997 data (4 points) and 
1999-2000 data (4 points). 

— The DVCS cross section differential in t, [56], Table 3a, (12 points). 

• Measurements by the ZEUS collaboration of 

— The DVCS cross section differential in t, [57], Table 1, only > 4GeV^ points (5 
points) . 

— The total DVCS cross section, [57], Table 4 (4 points). 

— The total DVCS cross section, [SS], Table 1 (6 points). 

In total, we have 95 data points. Fitted to these points was a version of the model used in 
[38] • Here we only describe differences, and list the free fitting parameters. The first difference 
is that, beside leading partial wave in S0(3) expansion of conformal moments of GPDs H^'^'^ 
and H*^ (the normalization of which is fixed by DIS F2 data, and residual t-dependence of H^'^^ 
being determined by free dipole mass parameter M^^^) and subleading partial wave with relative 
strength parameters s^'"^ and s^, we take here into account also the third partial wave with two 
new strength parameters s^^^ and sf. The second difference to [3S] is that ^tS is here modelled 
by a shape suggested by the pion pole contribution, but with normalization r^r and additional 
t-slope dipole mass as two additional free parameters. Together with parameters 
and 6™' parametrizing Jm'H, C and Mc parametrizing subtraction constanto, and f™'. 



^Note that the subtraction constant here is given by -2?, appearing in (|110p . 
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Table 3: Valence (top) and sea quark (bottom) related hybrid model parameters, extracted from the 
global DVCS fit KMM12. 

and parametrizing Um'H in the same way as in [3H], this brings the number of parameters to 
a total of 15. The fit of this model to all of the above data results in l^d.o.i. = 124.1/80, which 
is strictly speaking not a good fit, but it is acceptable for a global fit to data coming from such 
a variety of experiments and observables. Parameters of the fitted model are given in Tab. [31 
and the values resulting from this global fit for the 10 Hermes observables used for local fits in 
previous sections are given by solid lines on Figs. I5HTT] in Appendix \^ 

We now finally considerthe possible tension between our simple GPD model fit and the 
data. As one can see in Fig. |8l our fit reasonably describes the beam spin and beam charge 
asymmetries, but shows a slight tendency for the mean values to slightly overshoot the data 
values. This, as discussed above, is related to the problem of overshooting the beam spin 
asymmetry measurements, taken from events selected by the missing mass technique, with GPD 
models which employ a leading order description that have skewness ratios r > 1. The tension 
that appears in the lowest xb and bin for the first harmonic of the beam charge asymmetry 
may be related to a technicality; namely, our input scale for evolution is chosen to be 4 GeV^ 
and backwards evolution to < 1.5 GeV^ in the flavor singlet sector is a delicate procedure 
that is rather sensitive to the initial conditions. The dominant longitudinal proton spin-flip 
asymmetry is also overvalued by the model, see Fig. |9l which just reflects the fact that the 
description of unpolarized HALL A cross section measurements requires a large real part in 
the DVCS amplitude which is partially also obtained via a GPD dispersion relation from the 
imaginary part of CFF "H. Such a parametrization, which is effective for unpolarized DVCS 
measurements, is clearly disfavored if longitudinal proton spin-flip asymmetry data are included. 
The transverse proton spin- flip asymmetries, shown in Fig. [10] are well described by our Jm £ = 
3xnS = models, within the large experimental uncertainties. Mostly, this is also the case for 
the transverse proton spin-flip asymmetries, presented in Fig. [11] As an exception on might 
view given in f[36|) . where one may see some hint that the sign of the real part of £, 

adopted from the pion pole, contradicts data. Let us emphasize that the basic modeling of S in 
terms of pion pole contribution is oversimplified and already contradicts the GPD interpretation 
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of Tr"*" electroproduction data in the coUinear framework [l3l [59] . 

Compared to previous good DVCS world data fits to data from an unpolarized proton, having 
X^/d.o.f. ~ 1, we may interpret our findings here as a shght tension between a very simple 
model and fixed target data. One may argue that this tension is induced by the attempt to 
explain the unpolarized photon electroproduction cross section data from JLab Hall A by means 
of the common four twist-two GPDs, which have a oversimplified functional form. A definite 
conclusion cannot be given, however, since one may use more intricate parameterizations of these 
four twist-two GPDs, or use an improved framework by inclusion of perturbative corrections, 
certain twist-three [601 EB [62] and twist-four contributions [63| [M] , or extend the number of 
GPDs and thus CFFs. Certainly, there are also experimental uncertainties related to the issue 
of exclusivity, in particular the potential inclusion of events including a A-resonance within the 
experimental data set. We emphasize, however, that the inclusion of polarized proton data does 
not contradict the conclusion that GPD H plays the dominant role in the description of present 
DVCS data. 

4 Summary and outlook 

In this article we analyzed, by means of mapping and regression methods, the final set of DVCS 
off-t he-proton data from the Hermes collaboration extracted using a missing-mass event selec- 
tion method. Thereby, we still utilized the twist-two dominance hypothesis and, thus, we re- 
stricted ourselves to an overcomplete set of fourteen asymmetries. We showed that the Hermes 
collaboration provided an experimental proof of principle that, with an (almost) complete mea- 
surement of eight first harmonic asymmetries in the charge-odd sector, all four twist-two asso- 
ciated CFFs can be accessed in a BH dominated regime. In the remaining six zero harmonic 
asymmetries the CFFs are kinematically suppressed and we used them to check the twist-two 
dominance hypothesis or, alternatively, employed them to access twist-two associated CFFs. 
Higher harmonics, in the first place related by the remaining set of eight CFFs, were not con- 
sidered. 

Apart from restricting ourselves to the twist-two sector, our analysis is rather general. We 
tested different methods to access CFFs: mapping, local regression analysis, and model fits. In 
any CFF (or GPD model) extraction procedure from a set of asymmetries, one must bear in mind 
that separate solutions may exist for the BH dominated and the DVCS dominated regimes. Based 
on experimental evidence and model expectations, we took in our analysis the BH dominated 
solution. A one-to-one map of asymmetries to the space of CFFs can be either found by picking 
up the appropriate roots of a non-linear equation system (in our case eight equations that can 
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Figure 7: Left: Two projections witli total errors 0.01 ^ 6A < 0.03 (filled triangles) and 6A ^ 0.01 
(empty triangles) of charge odd asymmetries (left) for = 0.1, t = — O.lGeV^, and = 2.5 GeV^ 
are compared with Hermes measurements in bin ^2 (solid circles). Right: the resulting CFFs from a 
one-to-one map, based on the twist-two dominance hypothesis. 

be linearized) or, equivalently, by the least squares method. However, numerical noise in a blind 
"fitting" technique can provide an unphysical solution. Constraining the value of CFFs to force 
a "black box" fitting routine to provide the physical solution, a rather popular method, may 
have the disadvantage that a true one-to-one map cannot be obtained. The inclusion of more 
(twist-two related and unrelated) asymmetries in a fitting procedure will increase the noise and 
the physical solution is not necessarily obtained from the global minimum. 

In a one-to-one map of interference dominated, mostly charge-odd asymmetries, it turns out 
that only the imaginary part of CFF "H is not compatible with zero, and all other seven sub- 
CFFs can be considered as compatible with zero. Surprisingly, it turns also out that the best 
constrained quantity is the imaginary part of CFF Ti, followed by the real part and imaginary 
parts of "H, and to a much lesser extent also by the imaginary part of CFF S. The remaining four 
sub-CFFs are very poorly constrained. One may be tempted to consider the extraction of CFFs 
as a regression problem, and thus use selection criteria for filtering out the noise. As expected, in 
this approach it turns out that Jm "H is a robust quantity, the method does not necessarily yield 
a unique solution. We present two solutions: one obtained by the strict use of selection criteria 
and the other a more hand-picked solution informed by standard GPD model considerations. 
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The first one suggests that the negative real part of CFF £ can be considered as large and 
the second solution suggests that one may consider the real part of CFF "H as negative and the 
imaginary part of "H as positive. Certainly, from the GPD model point of view one would give the 
hand-picked standard solution preference. However, we emphasize that our one-to-one map only 
shows that this well constrained quantity JmH is small and does not allow a further partonic 
interpretation, e.g., extraction of any t-dependence. A rather analogous situation appears for the 
imaginary part of CFF £ — its magnitude is only loosely constrained by the data. The relation 
of GPD E to the spin sum rule inspires both experimentalists and theoreticians to make definite 
GPD model statements that certainly cannot, at present, be justified by simply describing 
experimental DVCS data with a given model. This we have illuminated in our global GPD 
model fit example where the GPD E (in the standard double distribution representation) has 
been set to zero, giving us the perhaps best (however, not perfect) DVCS world data description 
that is presently available. 

Let us illustrate that with the Hermes experiment, switched off 2007, our knowledge about 
CFFs could have been much better. Supposing that the longitudinally polarized proton asym- 
metries AuL,i and would have been measured in the charge-odd sector with the same 
statistics and slightly better systematics than the existing data taken with positron beam, the 
proton helicity conserved CFFs "H and "H could have been accessed with roughly the same accu- 
racy. Alternatively, a decrease of the large statistical uncertainties for the Alt,i quantities to a 
typical value of other Hermes asymmetries would have drastically improved the constraints for 
the real parts of proton helicity flip CFFs £ and £ (note, however, that these CFFs naturally 
suffer from a larger uncertainty). The effect of a decrease in the uncertainties that would be 
feasible at a fixed target experiment such as Hermes is illustrated by the filled triangles in 
Fig. [71 where we naively assumed for all polarized charge odd asymmetries a 5A ^ 0.03 error 
(statistical and systematic uncertainties are added in quadrature). 

The Hermes collaboration provided a total set of thirty-four asymmetries that includes an 
(almost) complete measurement of the second harmonic asymmetries in the charge odd sec- 
tor, which are primarily associated with four twist-three and four twist-two related CFFs. In 
principle, twenty-four a,sym.m.eiiies can be utilized in a one-to-one map, which leaves us ten asym- 
metries for a consistency check. Due to kinematic suppressions and noise it is expected that 
access to the eight twist-tree and transversity associated CFFs is not achievable; however one 
may hope for a test of whether our (and other) results are robust in an unbiased map, i.e. where 
one drops the twist-two dominance hypothesis. Before one undertakes such an attempt, several 
improvements and technicalities should be taken into account, which are partially also needed 
on more general grounds: 
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• A physically motivated parametrization of the DVCS tensor in terms of twelve CFFs. 

• A code which relates the set of twelve well defined CFFs to observables in an exact manner, 
perhaps for different conventions. 

• The (small) differences between proton polarization vector, defined in the lab frame, and 
those in the rotated frame, see Fig. [H should be taken into account, at least from the 
principled point of view. 

• It would be desirable that the covariance matrix for experimental measurements would be 
available for each given kinematical bin. 

Let us add to the first two points that various codes are presently used that relate (twist-two 
associated) CFFs (in most cases understood as quantities that are expressed by conventionally 
defined GPDs) that are based on different conventions and approximations. At present numerical 
differences do not matter, nevertheless, it is not desirable to proceed in this manner to the next 
generation of DVCS experiments. The remaining obstacles can be overcome by means of the 
parameterization for the virtual Compton scattering tensor and the analytic cross section results 
that are presented in [17]. 

In planned DVCS experiments at JLAB@12GeV, having an electron beam and fixed proton 
target, it is expected that experimental uncertainties will become much smaller, as already 
demonstrated by the hall a collaboration in measurements of DVCS cross sections. Altogether 
in such an experimental set up one can measure eight (seven) azimuthal angular dependent cross 
section combinations (asymmetries), compared to sixteen (fifteen) for both charges of electrons. 
Having precise cross section measurements, one can form, at least in principle, 4x4 even and 
4x3 odd harmonics, giving us twenty-four independent and four dependent (constant term) 
observables. Hence, one may have a handle on the CFF separation via a harmonic analysis |14j . 
Qualitatively or semi-quantitatively, one may also employ Rosenbluth separation to address 
the three separate parts of the photon electroproduction cross section. To which extent this 
separation is feasible requires detailed studies, where, unfortunately, the results will depend on 
assumptions. However, the possibility of obtaining any unbiased extraction of CFFs cannot be 
clearly stated at present. Certainly, a clean separation of the charge odd (interference term) and 
charge even part (sum of BH and DVCS cross sections) requires a high- luminosity positron beam, 
too. Supposing that in a (next-to-) next generation of lepton-proton scattering experiments both 
kinds of electrons are available and the total experimental uncertainties for all asymmetry (or, 
better, cross section difference) measurements are three times smaller than those of the beam- 
spin asymmetry measurements at Hermes, see empty triangles on Fig. [TJ we may expect that 
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the CFF 8 can be accessed in the BH-dominated DVCS regime with an accuracy that is presently 
available for CFF "H. 
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A Visualization of HERMES data descriptions 




Figure 8: Fits to liarmonics of asymmetries of scattering on an unpolarized target. Black dots are 
Hermes data with systematic errors added in quadrature. Local fits in two different scenarios are 
shown as red diamonds (fit to 3m7i and yitS) and blue pluses (fit to JmH, 9\eT-l, and 3mT-l), slightly 
displaced to the right for legibility. For comparison, we also show the result of a global fit to world 
DVCS data as a green solid line. 
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Figure 9: Fits to longitudinally polarized target asymmetry harmonics. Legend is as for Fig. [8] 
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Figure 10: Fits to transversely polarized target asymmetry iiarmonics. Legend is as for Fig. [8] 
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Figure 11: Fits to asymmetry harmonics of polarized electron/positron beam scattering on a transver- 
sally polarized target. Legend is as for Fig. [5] 
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